Model analysis of the world data on the pion transition form factor 
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We discuss the impact of recent Belle data on our description of the pion transition form factor 
based on the assumption that a perturbative formalism and a nonperturbative one can be matched 
in a physically acceptable manner at a certain hadronic scale Qo- We discuss the implications of the 
different parameters of the model in comparing with world data and conclude that within experi- 
mental errors our description remains valid. Thus we can assert that the low Q 2 nonperturbative 
description together with an additional 1/Q 2 term at the matching scale have a strong influence on 
the Q 2 behavior up to very high values of Q 2 . 

PACS numbers: 12.38.Lg, 12.39.St, 13.40.Gp, 13.60.Le 

New data of the pion transition form factor (ttTFF) from the Belle collaboration have just appeared 
These data, above 10 GeV 2 , are smaller in magnitude than the previous BABAR data [J, which 
generated considerable excitement. The question to unveil is the scale of asymptotia. BABAR data, taken 
at face value, implied that asymptotic QCD behavior lies at much higher Q 2 than initially expected 
Belle data seem to lower that scale. We show here that our scheme can accomodate easily all data without 
changing the physical input. 

At the time of the BABAR data we developed a formalism to calculate the irTFF Q , which consists of 
three ingredients: i) a low energy description of the nTFF; ii) a high energy description of the irTFF; 
Hi) a matching condition between the two descriptions at a scale Qo characterizing the separation between 
the two regimes. For the low energy description we took a parametrization of the low energy data to avoid 
model dependence at Qo- The high energy description of the irTFF, defined by the pion Distribution 
Amplitude (ttDA), contains Quantum Chromodynamic (QCD) evolution from Qo to any higher Q, a 
mass cut-off to make the formalism finite, and an additional 1/Q 2 term which leads to modifications of 
the matching condition. 

Let us recall some aspects of the formalism. The high energy description, to lowest order in perturbative 
QCD, for the transition form factor in the process ir° — > 7 7* in terms of the pion distribution amplitude 
(ir DA), is given by 

Q 2 F(Q 2 ) = ^f -%tM*, Q 2 )- (i) 

a Jo x + -qt 

We follow the proposal of Polyakov Q and Radyushkin Q and introduce a cutoff mass M to make the 
expression finite. Q 2 = —q 2 , is the momentum of the virtual photon, <p v (x, Q 2 ) is ttDA at the Q 2 
scale and f v = 0.131 GeV. In this expression, the Q 2 dependence appears through the QCD evolution of 
the it DA. 

Despite the fact that several models reproduce the low energy data, in order to have a model inde- 
pendent expression for the form factor at low virtualities, we adopted a monopole parametrization of the 
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TtTFF in the low energy region as 
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FIG. 1: We show the result for the transition form factor in our formalism for M = 0.690 GeV, a — 0.032 and 
Cz = 2.98 10~ 2 GeV 3 and defining the matching point at Qo = 1 GeV (solid line). The band region results from 
the indeterminacy in Aa — ±0.004. The lower plot shows the detailed behavior for low virtuality. Data are taken 
from CELLO [Tof, CLEO [H, BABAR [| and Belle Q]. 



Additional power corrections can be introduced in Eq. [T] by adding to the lowest order calculation a 



3 



0.3 

o 
O 

•C 2 
0.1 



10 20 30 40 

Q 2 (GeV 2 ) 




2 4 6 8 10 12 



Q 2 (GeV 2 ) 




FIG. 2: We show the result for the transition form factor in our formalism for M = 0.620 GeV, a — 0.032 and 
the value of C3 = 1.98 10~ 2 GeV 3 corresponding to 20% of the contribution at the matching point at Qo = 1 
GeV (solid line). The band region gives the variation of the results due in ±10% in the contribution of higher 
twist. The lower plot shows the detailed behavior for low virtuality. Data are taken from CELLO [1Q], CLEO 
[H, BABAR and Belle Q]. 



term proportional to Q 2 , 

Q F iQ ) = —5— / Tf2^{x,Q ) + — j . 

6 Jo x + |p- Q 

Using a constant it DA the matching condition becomes [f|, 
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with Qo = 1 GeV. This equation allows to determine M, once we have fixed the value of C3. 
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FIG. 3: We show the result for the transition form factor in our formalism for M = 0.690 GeV, a = 0.032 and the 
value of Cz — 2.98 10 -2 GeV 3 corresponding to 30% of the contribution at the matching point at Qo = 1 GeV 
(solid line) .The lower plot shows the detailed behavior for low virtuality.The dotted curve represents the higher 
twist contribution.Data are taken from CELLO [3, CLEO BABAR @] and Belle 



We analyze here the sensitivity of the data to the various parameters involved. We keep as close 
as possible to our previous fit analyzing the data with respect to small variations in the low virtuality 
parameter a and in the higher twist parameter C3. In FigfJJwe show the effect of the precision in the 
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determination of the monopole parametrization. We see that as a increases from 0.032 to 0.036, i.e. 
within the error bars, the irTFF decreases. The sensitivity to C3 is shown in Fig. [2] and we note that 
as the value of C3 increases from C3 = 0.99 10 _2 GeV 3 , which corresponds to a 10% contribution to the 
form factor at Qo, to 2.98 10 -2 GeV 3 , which corresponds to a 30% contribution, again the value of the 
TtTFF decreases. Thus a small increase in a and C3 moves our result toward the Belle data. Finally, in 
Fig. [3] we plot the better fit (x 2 /dof — 1.21) taking into account all the world data which corresponds 
to a = 0.032 with the C3 term at the 30% value. We stress that there is no strong correlation between a 
and C3 as long as a is kept within its experimental error bars. Thus the fit is quite stable with respect 
to the parameters of the low energy model. 

The fit to the data is excellent with a very small variation of the 1 /Q 2 contribution at Qq from previous 
fit, i.e. from 20% to 30%. It must be said, before entering the discussion of this fit, that in our previous 
work 5j] we pointed out that the average value of the hig hest energy data points of BABAR were too 
large, a conclusion reached also by other analyses fi~2l. Il3|. In Fig. [3] we show not only the fit for 30% 
contribution of C3/Q 2 at Q , but its behavior for higher values of Q 2 . As can be seen, also stressed in our 
previous work, this contribution is small in size. However, and this an important outcome of our analysis, 
it is instrumental in fixing the initial slope at the matching point, which determines, after evolution, the 
high energy behavior of the form factor. 

In our opinion the Belle data confirm the BABAR result that the irTFF crosses the asymptotic QCD 
limit. This limit is well founded under QCD assumptions, but nothing is known of how this limit is 
reached, if from above or from below. BABAR and Belle data suggest that the limit is exceeded around 
10 — 15 GeV. Our calculation is consistent with this result. The necessary growth of the irTFF between 
5 — 10 GeV to achieve this crossing is in our case an indication of nonperturbative behavior and C3/Q 2 
contribution at low virtuality. The determination of the crossing point is a challenge for any theoretical 
model and therefore, the precise experimental determination of it is of relevance. Many models fail to 
achieve this crossing because their pion DA is defined close to its asymptotic form. 

The pion DA can be expressed as a series in the Gegenbauer polynomials, 

<t>, (x, Q 2 ) = 6x (1 - x) 1 1 + J2 a n (Q 2 ) C Z J 2 (2x - 1) J (5) 

\ n(euen)- 2 / 

We can compare different models by looking at the values of the coefficients of the expansion a n (Q 2 ) ■ 
In our case, at Q 2 = 1 GeV 2 many a n coefficients are significant, but we focus our attention in a few 
terms: a 2 = 0.389, a 4 = 0.244 and a 6 = 0.179. At Q 2 = 4 GeV 2 we obtain the values a 2 = 0.307, 
a 4 = 0.173 and ag = 0.118, which are close to those obtained by Polyakov |(|. Consistently, our result for 
the irTFF is similar to that obtained in ref. [6]. At Q 2 = 5.76 GeV 2 we obtain a 2 = 0.292, a 4 = 0.161 
and a6 = 0.108, which are very different from those of ref. These author use for their fit BABAR 

data for the ijTFF [Hj], together with the pion data. It is therefore not a surprise that these authors come 
to a different conclusion, namely, that the Belle and the BABAR data cannot be reproduced to the same 
level of accuracy within the Light Cone Sum Rules approach [l6| . However, in an extension of the ideas 
developed in the present paper to the r\ case studied in ref. [17] looking at the state \q) — | (\uu) + \dc£)) 
a very different structure of the a n coefficients to that of the pion arises. At Q 2 — 1 GeV 2 , the values of 
the coefficients are a 2 = 0.134 and a 4 = 0.352 or, equivalently, at Q 2 = 5.76 GeV 2 we have a 2 = 0.101 
and a 4 = 0.232. Therefore, that study does not supports the combined use of both data sets 

We have developed a formalism to describe the irTFF on all experimentally accessible range, and 
hopefully beyond. The formalism is based on a two energy scale description. The formulation in the low 
energy scale is nonperturbative, while that of the high energy scale is based o n p erturbative QCD. The 
two descriptions are matched at an energy scale Qo called hadronic scale [IH, [lj| . We stress the crucial 
role played by the nonperturbative input at the level of the low energy description. It is an important 
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outcome of this calculation the role played by the 1/Q 2 power correction term in determining the slope 
of the data at high Q 2 , despite the fact that they do almost not contribute to the value of the irTFF . 

We have used a flat n DA, i.e. a constant value for all x @, 0] , which with our normalization becomes 
4>{x) = 1. Our choice has been motivated by chiral symmetry [f|. Model calculations, Nambu-Jona- 
Lasinio (NJL) [20l - 423j and the "spectral" quark model give a constant tt DA. The ttTFF calculated 
in these models, however, overshoots the data [25| . emphasizing the importance of QCD evolution. 

The calculation shown proves that the BABAR and Belle results can be accommodated in our scheme, 
which only uses standard QCD ingredients and low energy data. Moreover, at the light of our results, 
we confirm that at 40 GeV 2 we have not yet reached the asymptotic regime which will happen at higher 
energies. 
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